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Introduction 
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This is the final report which summarizes the research accomplishments under the 
project entitled "Physics of Heat Pipe Rewetting" under NASA Grant No. NAG 9-525, 
Basic, during the period of April 1, 1991 to January 31, 1994. 

The objective of the research project was to investigate both analytically and 
experimentally the rewetting characteristics of the heated, grooved plate. The grooved 
plate is to simulate the inner surface of the vapor channel in monogroove heat pipes for 
space station design. In such designs, the inner surface of the vapor channel is threaded 
with monogrooves. When the heat pipe is thermally overloaded, dryout of the 
monogroove surface occurs. Such a dryout surface should be promptly rewetted to 
prevent the failure of the heat pipe operation in the thermal radiator of the space station. 

Experimental Setup 

The schematic of the experimental setup is shown in Fig. 1. The test section is a 
grooved, oxygen-free copper plate whose dimensions and design are shown in Fig. 2. 
The plate temperature was monitored by 16 copper-constantan thermocouples mounted 
at 1.2 mm below the upper surface of the plate (see Fig. 2b). While these thermocouples 
provided plate temperatures at selected locations, a traversing surface measurement 
device was mounted above the test section to provide a continuous measurement of the 
surface temperature axially along the plate. The plate is heated and maintained at a 
constant temperature at one end by a heater assembly. Then the coolant, isopropanol 
liquid, is pumped to the reservoir and to wet the other end of the hot grooved plate. 
The plate is quenched by the advancing isopropanol liquid film in the grooves. The 
transient axial temperature profiles of the plate were recorded by the data acquisition 
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system. The wet front velocity was recorded by a personally owned 8 mm video 
camcorder. However, when the tape of the camcorder was played back frame by frame, 
it was difficult to determine the wet front position versus time since the picture was 
blurred with streak lines. A better, high precision camcorder and a compatible four- 
head VCR are apparently needed to determine the wet front velocity. Due to the 
limited funding provided, no budget was available for such purchases. 

Experimental and Theoretical Results 

Rewetting experiments of the heated, grooved plate were performed with the plate 
placed horizontally and then placed with an upward inclined angle to see gravity effects. 
Figure 3 shows the comparison of the temperature profiles of the horizontal plate and 
the plate with an upward inclined angle of 30°. The temperature profiles of the inclined 
plate drops slower than the horizontal plate, indicating that the gravity effect is to slow 
down the advance of the rewetting front. Figure 4 shows the comparison of the 
temperature profiles at various times of the horizontal plate which was placed in face 
up (namely, the grooved side of the plate was facing upward) and in face down 
positions. Again, the difference in the temperature profiles shows that gravity has an 
effect on th6 rewetting of the grooved plate. 

Theoretical investigation of the rewetting phenomena has also been conducted. The 
results of the study have been published in The Journal of Heat Transfer. A copy of the 
reprint is attached at the end of this report for completeness. In the paper, a model was 
presented to analyze the rewetting process of the plate with and without grooves. 
Approximate analytical solutions were presented for the prediction of the rewetting 
velocity and the transient temperature profiles of the plate subject to uniform heating. 
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Also presented is a simple method to predict the dryout heat flux of a liquid film 
flowing over a heated smooth and grooved plate. The theory presented in the paper is 
applied to predict the transient temperature profiles of the present experiments. Since 
the heated plate is not subject to a uniform heating but rather, the plate is heated at one 
end, no analytical solution seems obtainable using the model proposed in the paper. 
Instead, numerical solutions using the model presented in the paper were obtained and 
compared with the experimental data as shown in Figs. 5 and 6 for the rewetting of the 
hot grooved plate initially at 150°C after t = 20 and 40 seconds, respectively. The 
comparison shows that the actual temperature in the dry region (i.e., the region ahead 
of the rewetting front) is lower than the prediction. This indicates that there was heat 
loss in the dry region that was not accounted for in the model. On the other hand, the 
actual temperature profile in the wetted region drops slower than the prediction. This 
is possibly caused by the two-dimensional conduction effects in the plate which are not 
accounted for in the model. Nevertheless, the major feature of the rewetting 
temperature profile is satisfactorily predicted by the model. 
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(a) Grooved Microfin Surface 


(b) Thermocouple Probe in Plate 



(c) The End of Grooved Plate 


Fig. 2. 


Schematic of the grooved copper plate. 




COMPARISON BETWEEN THE RUN UP AND RUN DOWN TEST 

(Gravity Testing, T(1]=300 C) 
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plate in face up and face down conditions. 



: Experimental data 
: Numerical results 
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Comparison of predicted and measured temperature profiles of the grooved 
plate heated one-end at 150°C. The profile is 20 seconds after rewetting. 



160.00 


UD 



CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

lO 

CD 

CD 

HD 

CD 

QnI 

CD 

j 

CD 

CD 

CD 

CD 

CD 

CO 


r- j k ^ — 

07 3jn^DJ9duJ3i 


Comparison of predicted and measured temperature profiles of the grooved 
plate heated one-end at 150°C. The profile is 40 seconds after rewetting. 
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Rewetting Theory and the Dryout 
Heat Flux of Smooth and Grooved 
Plates With a Uniform Heating 

The evaporation and condensation of thin liquid films are of significant importance 
in a wide variety of problems ranging from specific applications in the heat pipe 
field to more general ones in the chemical t nuclear , and petrochemical industries. 
Although several investigations have been conducted to determine the rewetting 
characteristics of liquid films on heated rods , tubes, and flat plates , no solutions 
are yet available to describe the rewetting process of a hot plate subjected to a 
uniform heating. A model is presented to analyze the re wet ting process of such 
plates with and without grooves. Approximate analytical solutions are presented for 
the prediction of the rewetting velocity and the transient temperature profiles of the 
plates. It is shown that the present rewetting velocity solution reduces correctly to 
the existing solution for the rewetting of an initially hot isothermal plate without 
heating from beneath the plate. Numerical solutions have also been obtained to 
validate the analytical solutions. Finally , a simple method is presented to predict 
the dryout heat fux of a liquid film f owing over a heated smooth or grooved plate. 
The results of the prediction are found to be in reasonable agreement with the 
existing experimental data. 


Introduction 

The rewetting process is a conjugated heat transfer problem 
involving interactions between a solid wall and flowing fluids. 
The process for rewetting of a grooved plate with a uniform 
heating is complicated, as the rewetting velocity varies with 
time, physical geometry of the grooves, plate properties, fluid 
properties, and the applied heat flux. Although several inves- 
tigations (Yamanouchi, 1968; Thompson, 1972; Duffey and 
Porthouse, 1973; Sun et ah, 1974; Alario et al., 1983; Grimley 
et al., 1988; Stroes et al., 1990; Ferng et al., 1991; Peng and 
Peterson, 1991) have been made to determine the rewetting 
characteristics of liquid films on heated rods, tubes, and flat 
plates, none has yet presented the solution for the rewetting 
process of a heated plate with a smooth or grooved surface 
subjected to a uniform heat flux. Is is noted that the surface 
with small grooves has received increasing attention as it has 
many practical applications. For instance, the microgrooved 
surface is employed most often to enhance heat transfer (Grim- 
ley et al., 1988). Microgrooves are also useful for replacing 
the wicking material in heat pipes. In fact, they are used in 
the innovative monogroove heat pipe design for the thermal 
radiators of the space station (Alario et al., 1983). When the 
thermal radiator is overloaded with a heat flux discharged from 
the condenser of the thermal bus system in the space station, 
the circumferential section of the monogrooved heat pipe di- 
rectly underneath the heat flux is the location where dryout 
occurs first. In a recent space shuttle flight test, dryout of the 
heat pipe did occur. Therefore, it is of interest here to inves- 
tigate the rate of rewetting speed and the condition that leads 
to the dryout of a plate subjected to uniform heating. 

There are fundamental differences between the rewetting of 
a hot plate with a uniform heating and that without heating. 
As shown later (Eq. (1) and Appendix A), unlike the case 
without heating, the heat conduction equation is a transient 
equation even after the equation is transformed to the La- 
grangian coordinate moving with the liquid rewetting front. 


Contributed by the Heat Transfer Division and based on a paper presented 
at ASME Winter Annual Meeting, Anaheim, CA, November 8-13, 1992. Man- 
uscript received by the Heat Transfer Division October 1992; revision received 
June 1993. Keywords: Augmentation and Enhancement, Heat Pipes and Ther- 
mosyphons, Space Power Systems. Associate Technical Editor: L. C. Witte. 


Secondly, the plate temperature in the dry, insulated region 
far ahead of the rewetting front is not only invariant in the 
axial direction, but also a linear function of time due to a 
uniform heating (see Eq. (8)). A recent attempt (Peng and 
Peterson, 1991, 1992) to solve for the rewetting process of a 
plate subjected to a uniform flux has not accounted for these 
differences. They apparently encountered a compatibility dif- 
ficulty in matching the solution of the temperature profile of 
the heated plate with a proper boundary condition because of 
the use of the steady-state transformed heat conduction equa- 
tion, which is basically valid only for the case without internal 
or external heating (Yamanouchi, 1968; Thompson, 1972; 
Duffy and Porthouse, 1973; Sun et al., 1974; Tien and Yao, 
1975). 

One of the objectives here is to present a physical model 
suitable for the rewetting analysis of a plate heated by a con- 
stant heat flux from below. Solutions for the transient tem- 
perature profile and the rewetting velocity are presented. The 
incompatibility problem is also resolved. Another objective is 
to analyze the dryout limit of a grooved surface initially wetted 
by a liquid film. A means is provided for the prediction of the 
maximum heat input that results in the dryout of the liquid 
on the plate. The result of the prediction is compared with the 
experimental data of Grimley et al. (1988). 

Rewetting Model and Solutions 

We consider the rewetting process of a hot plate initially at 
a uniform temperature T { with no liquid on the plate as shown 
in Fig. 1. The plate is heated from below by a uniform heat 
flux and is quenched by a liquid advancing along the direction 
of the grooves on the top surface of the plate. In order to 
simplify the complexity of the physical phenomena of the re- 
wetting process, we consider first a heated smooth surface plate 
with no parallel grooves. 

1 Smooth Surface Plate. The rewetting process of a hot 
dry plate is sketched in Fig. 1. A liquid film from a liquid 
reservoir, driven by its surface tension, is to advance along the 
hot plate. Similar to prior studies by Yamanouchi ( 1968 ), Duf- 
fey and Porthouse ( 1973 ), and Sun et al. ( 1974 ), the initial 
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temperature of the plate is assumed to be higher than the 
Leiden frost temperature T 0 such that the rewetting process is 
assumed to be conduction controlled. For a liquid film on a 
plate, there exist three heat transfer regimes: boiling, convec- 
tive evaporation, and single-phase convection. However, sim- 
ilar to these prior studies, the following assumptions are made: 
a constant averaged heat transfer coefficient, for simplicity, 
in the wet region to remove the heat from the thin plate to the 
liquid film, no heat loss to the environment in the dry region, 
the plate at the rewet front remains at a constant Leidenfrost 
temperature, the liquid film remains at the saturation tem- 
perature T sy and the plate is thin enough that the one-dimen- 
sional rewetting model can be invoked. It is therefore proposed 
to solve the following governing equation (see Appendix A) 
on a Lagrangian coordinate moving with the rewetting front: 


where 


P= 


A = 


so d 2 e 

do 
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where B is the Biot number; A y P t 6( 77 , r), 77 , and r are the 
dimensionless heat source, rewetting velocity, temperature, 
length, and time respectively. In the wet region (~L X < x < 
0 ), h - const ^0 while in the dry region (0<x<Z, 2 )» h = B-0. 

The above governing equation is different from the equation 
in all prior rewetting models (Yamanouchi, 1968; Du f fey and 
Porthouse, 1973; Sun et al., 1974; Peng and Peterson, 1991, 
1992) in that an extra term, namely, the transient term on the 
left-hand side of the equation, is added. The addition of this 



Fig. 1 Schematic of the rewetting film on the heated p| a t e 


transient term is essential in analyzing the rewetting proccs\ 
of the heated plate subjected to a heat flux conduction asshoun 
in Fig. 1 . This is because the temperature profile of the piar c 
in the Lagrangian coordinate (x, y) is no longer invariant with 
time as in the case of rewetting an initially hot isothermal pb rc 
without heating from beneath the plate. Without it, the in 
compatibility difficulty as noted above will arise and the Tinal 
solution cannot be expected to satisfy the transient boundary 
condition at x= 00. The addition of the transient term is there 
fore essential. It, however, renders some mathematical com 
plications in the solution of the rewetting velocity. 

The initial condition is 


0(7, = 

y 0~ *s 

while the boundary conditions are 


t) = 0 

0(0, r) = 1 

e (lL 2 > 

J o~ * $ 

where i) U and y L2 are defined as 



(}) 


Hi 

(5) 

(6) 


(7) 


Nomenclature 


A = dimensionless heat source 
a\ = coefficient defined in Eq. 
(65) 

B = Biot number with respect to 
the convective heat transfer 
coefficient 

C p = thermal capacitance, J/kg- 
°C 

h = surface convective heat 

transfer coefficient, W/m 2 - 
°C 

h - average boiling heat transfer 
coefficient, W/m 2 -°C 
h a - convective coefficient be- 
tween plate and environmen- 
tal gas, W/m 2 -°C 
K = thermal conductivity, 
W/m-°C 

L = length of the plate, m 
L\ - length of wet region, m 
Li = length of dry region, m 
N = grooved geometric coeffi- 
cient 

n - exponent defined in Eq. (65) 
P = dimensionless rewetting 
velocity 


q = uniform heat flux, W/m 2 
<7chf = critical heat flux (CHF), 
W/m 2 

Qi = incipient boiling heat flux, 
W/m 2 

<7max = maximum heat flux, W/m 2 
s = plate thickness for smooth 
(5,) or grooved (s x - f,) plate, 
m 

s 1 = plate thickness, m 
s 2 = liquid film thickness, m 
T - temperature, *C 
T a — environmental temperature, 

•c 

T f - inlet liquid temperature, °C 
Tf = incipient boiling tempera- 
ture, °C 

To = Leidenfrost temperature, °C 
T s = saturation temperature, °C 
T, = initial hot surface tempera- 
ture, *C 
t = time, s 

U r = rewetting front velocity, m/s 
x = length of liquid film in mov- 
ing front coordinate, m 


x' = length of liquid film in sta- 
tionary system, m 
oq = constant defined in Eq. (23) 
a 2 — constant defined in Eq. (45) 

01 = constant defined in Eq. (23) 

02 = constant defined in Eq. (45) 
17 = dimensionless length coordi- 
nate with respect to x 

tj' = dimensionless length coordi- 
nate with respect to x' 

0 = dimensionless temperature 
B h — dimensionless temperature 
of the transient part 
8 S - dimensionless temperature 
of the steady-state part 
p = density, kg/m 3 
t = dimensionless time 

Subscripts 

d = dry region 
g = grooved plate 
L = total length of the plate 
L x = length of wet region 
L 2 = length of dry region 
s = smooth plate 
w = wet region 
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and T l is the dry plate temperature at 17^. The condition given 
in Eq. (6) implies that the plate is long enough that, within 
the time period of interest, the plate temperature at r ) U is not 
thermally affected by the rewetting process. Consequently, it 
can be expressed as 




<7 1 


( 8 ) 


pCpS x 

which is used in Eq. (6). 

In the wet region ^ 17 ^ 0), the mathematical model 

of this problem is given as 

nn 

(9) 


dd d z e n do _ 

dT~dr, 1 + P d V + 


e(i j, o) = e, 

0(-Vl,,t) = 0 


d~d s w dd s w 

— ^ + P — ^-Bd^ W + A =0 


dr) drj 

Os, w(0)=l (15) 

Os.U-vl^ 0 (16) 

while Of,, w is the solution to the transient part of the problem 


dVh, w d 2 6 hf w _ 30*, 




3r dq dr? 

0*. »(>), O) = 0|-0 J, „(t|) 

r) = 0 

0 *. J.~VL r T) = 0 

The solution of the steady-state problem is given as 

A_ 


(17) 

(18) 

(19) 

( 20 ) 


9s. J.v) = 


-F 


H) 




B 


e' r '”h -e' r ^h 


e' 


(-4 


1 r, ^i + — 


e - r Hz., _ e ~ r 2 , »e l 


e r * + - (21) 


n = 


r 2 = 


- P+\fp z + 


45 


— P—\J P 2 + 4B 


To solve the transient solution, the following transformation 
is introduced: 


r) 

where 

a, = -P/2, 0 l =-B-P 2 /4 
Then Eqs. (17) to (20) become 

3r 3i; 2 

»(?, 0) = e-^(^-^ w ( n )) S / l (i 7 ) 
u(0,t) = 0 
t) = 0 


( 22 ) 

(23) 

(24) 

(25) 

(26) 
(27) 


( 10 ) 

(11) 


The solution of the above can be readily obtained from Carslaw 
and Jaeger (1959) as follows: 


<9(0, r) = I (12) 

Solving for the exact, analytical solution of the above equa- 
tions appears to be difficult, if not impossible. Therefore an 
approximate, analytical solution is sought. This is made pos- 
sible by treating the Peclet number, P, as a constant value in 
the mathematical deliberation in order to achieve a closed- 
form solution. This approximation appears to be reasonable 
as the rewetting velocity tends to reach a quasi-steady state 
after an initial period when the liquid film is brought into 
contact with the hot plate. The numerically exact solution will 
be presented later to check the accuracy of the closed form 
solution. Accordingly, the solution is split into two parts, 

0(V, r) = 6 St M + Oh.M, t) (13) 

where is the solution for the steady-state part of the prob- 
lem, 

f2a . ... dQ. ... 

(14) 


where 


/»=! L ^IJ 

IJl, J- 


(28) 


Tim sin 


1 


MV)=e-°' v 0, 


1-4- 


J~ — «T V 
L 9l, 

H)*"" l,t 5 


dV (29) 


e-'iiLi-e-'m, 


Therefore, the combined solution is 
9(y, t) = 6 s „(v)+9/,, „(>?, t) 


B e’^x-e'^h 




*'”'-51 <30 > 


A 

’ b' 


H> 




e~ r ^t-e~ r ^h 




H) 


e~ r ' n L i+- 


" r l^l _ P ~ r VL, 






V- 1 ^ on 

n* 1 \ ^ L \) 


where a nl and/ t ( V) are given by Eqs. (29) and (30), respectively; 
a 1 and are given by Eq. (23). 

Before proceeding to the solution for the dry region, it is 
noted that the steady-state solution given above is identical to 
the solution presented by previous workers (Sun et al., 1974) 
when A -0 (or q = 0). 

Similarly, the problem in the dry region (0 < 17 < 1 j u ) is 
described by 


30 3 J 0 30 


d rt? +p ^ +A 

(32) 

0ft.O) = 0, 

(33) 

0(0, t)=1 

(34) 

0(iJc 2 . r) = 0, + At 

(35) 
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which is likewise split up into the steady-state and transient 
solutions. 




dOs. d 

dtf 1 dr) 

9,. A 0)=1 
9 S , d 0 Jt 2 ) = 0| 


+ /1 =0 


and 


d9h,a _ d 9h,_d + p 99jtjd 

dr dr? di) 

Oh.d (v. 0 ) = O I -0 S , M 
9,. AO, t) = o 

9h,d(VL 2 . t) = At 

The solution of the steady-state problem is 

A 
~P‘ 


(36) 

(37) 

(38) 


(39) 

(40) 

(41) 

(42) 


Os. <Xt) = 


t 'i + 77" 2 
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e~ mL i - 1 


A 

V 1 


„ A 

e ' + yit-i- 1 
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A similar transformation is introduced 

9>.d(v, T) 

where 

t*2 = -pf 2, &i = - p J /4 

to reduce Eqs. (39) to (42) to the form 
du d 2 u 

dr~ dn 2 

u(r,,O) = (0 l -8 S 'M)e- a »*M V ) 

w(0, t) — 0 

«(>j£ 2 . T) = /tT-e“ <a 2’i.2*42 r »S02(T) 


(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 


which can be solved by the method of linear superposition. 
Thus, the transient temperature profile in the dry region can 
be readily found as 


6(v, 7 ) = e“2’ + 4 2 rj’g^ sin e -" Vr/ ’t 2 

+ — (At) — V (- i)" Sm ^"7 

Vl 2 P„ 


Upon the substitution of Eqs. (31) and (50) into Eq. (51), the 
dimensionless rewetting velocity P can be determined by the 
expression 




+ 4B 






-'I 

=^«,2 2 
VW 1i. 2 

} 2 T 00 fylp- a 2*£.}-0*r . 


+ (0^-02)^ 


»‘ fl2,r -l) j 


. /4 


(52) 


It is of interest to examine the limiting solution of the above 
for the case that has been investigated by Yamanouchi (1968), 
namely, the rewetting of an infinitely long, hot, and isothermal 
plate without any heating. In this case, A =0 (i.e., g-0) and 
then by setting rj L] oo t r— oo, the above solution reduces 

to 




(53) 


r -*- 

J o 


tyd<f> 2 (X) + 1 


A 

( 9 > + p Vl 2 ~ 

— j 


(9\+~ i?c 2 -l^ 


-PVLj 


1 




(50) 


where 0„ = nr/i) L2 and 

2 f u 2 t mrw 

@n 2 ~ — \ / 2 (w) sin dw 

Vl 2 j o *7 l 2 

At the rewetting front, the conductive heat flux is continuous 
(Yamanouchi, 1968; Thompson, 1972; Duffey and Porthouse, 
1973; Sun et aL, 1974), i.e.. 


which is exactly the same as the well-known Yamanouchi’s 
solution (1968). 

Due to the absence of data on rewetting velocity on a hot 
plate heated by a uniform heat flux beneath the plate, the 
solution given by Eq. (52) cannot be compared with experi- 
mental data. However, experimental data are available on a 
hot plate without heating from below. 

Yamanouchi (1968), for example, has confirmed reasonable 
agreement between the limiting solution given by Eq. (53) and 
his data. The general solution presented above for a smooth 
plate will be used in the following section to yield the solution 
for the plate with axial grooves. 

2 Grooved Plate. The rewetting model of the grooved 
plate is based on that of the smooth surface plate. The coolant 
is driven by the wicking (surface tension) effect of microfins 
and is assumed, without loss of generality, to fill up the grooves 
as shown in Fig. 2. 

At the level of y — 0, for a thin plate and thin grooves, 

3T 

- K ^\,-o-(2? 2 + t 3 )~h(T-T s ) 

MKx+h) + h a (T-T a ).2! 2 (54) 

where h is the convective coefficient of a smooth surface, T a 
is the environmental temperature, h a is the convective coef- 
ficient between the plate surface and the environmental gas 
above (h a = 0). At the level of y= - (s, - f,). 


36 

dij 




31 , 

dy 




<7 

' K 


(55) 


(51) 


Since the plate is thin. 
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Fig. 2 Microfin structure and cross-sectional shape of grooved plate^ 


d 2 T 

— 2 = lim 

dy 


VdT _3T - 
dy y+iy dy y 


by 


3T 3T 
dy ym0 ~dy 




(56) 


Combining Eqs. (54), (55), and (56) yields 


d 2 T 

v 


K 


{T-Td-± 




(57) 


where the grooved geometric coefficient is defined as 

N=( 2f( + «/(% + « 

The result given by Eq. (57) suggests two useful simplifi- 
cations. First, the factor (Nh) is the equivalent convective heat 
transfer coefficient of the grooved plate and can be approxi- 
mated by that of a thin liquid film on a smooth plate multiplied 
by a factor of N, namely. 


/heat transfer coefficient 


of a grooved plate 




heat transer coefficient 1 ) 


of a smooth plate 


where N is more generally defined as 




the wetted perimeter 
width of the cross section 


(58) 


Second, the governing equation for the grooved plate with a 
uniform heating remains unchanged provided h is replaced by 
(Nh) and the dimensionless variables are properly scaled as 
follows: 

P = U,pC p (s x -t x )/K\ B = Nh(s x ~ r , )/K 
A —q(s x — ff) /[K( T 0 - T,)] ; rj = x/( s x - f,) 


Vli “ -f-i/ (-S' i — fi ); 


VL2 = W(S\ ~ f|) 


T=t/[(s i -e l ) 1 P c p /K] 


(59) 


where h is the convective coefficient of the smooth plate with 
a uniform heating. Therefore, the solution of the grooved plate 
with a uniform heating in the wet region is the same as Eq. 
(31) with the modification above and the solution for the dry 
region is identical to Eq. (50). 

Numerical solutions have also been obtained by solving for 
the original governing equation (see Appendix A, Eq. (A7)) 
fixed to the nonmoving coordinates (*', /), 



OmansJooi*$5 Length Cn> 

Fig. 3 Numerical solution for grooved plate 



Fig. 4 Rewetting velocity comparison, ^ = 2.5143, ^ = 14.5445, 
S = 0.0693 (wet region), A - 0.01954 (Case 1), A = 0.08142 (Case 2) 


subject to 


30 

3t 


d 2 e 


7-7 i-BO + A 


*(i?\ 0) = *, 

6 ( 0, t) = 0 
Q(vl> r) = 0, + At 


(60) 

(61) 

(62) 

(63) 


where rj' =x'/5, t)' l = L/s , and5 = Sj and -f,) for the smooth 
and grooved plates, respectively. The rewetting front location, 
Tj£ r is determined from 6 ( 77 ^. 1 » t) = 1. The differential equation 

is discretized by a standard finite-difference approach. Figure 
3 illustrates the computed wall temperature of the grooved (or 
smooth) plate versus the plate length, the so-called rewetting 
temperature curve, at various times when the liquid film is FC- 
72, which is Case 1 of Fig. 4. These numerical values are 
selected from the experimental condition of Fig. 5 of Grimley 
et al. (1988), namely, for FC-72 fluid on the smooth surface, 
7)~r,= irc, r 0 -7;=17.5 °C, 7>56°C, r, = 100°C from 
which Oi is calculated. Other dimensionless parameters A , B, 
and i} L are estimated from the following values. The needed 
heat transfer coefficient h — 2779.49 W/m 2 C is calculated from 
their boiling curve by Eq. (65). For the grooved copper plate, 
C p = 383.1 J/kg°C, p = 8954 kg/m 3 , K = 386 W/m*C. When 
the liquid fills up the grooves as shown in Fig. 2, the grooved 
geometric coefficient is N= 1.75, which is based on their geo- 
metric dimension of (\ = 0.5 mm, t 2 = 0.2 mm, f 3 = 0.4 mm, s x = 6 
mm, and q = 24,000 and 100,000 W/m 2 for cases 1 and 2, 
respectively. They also reported dryout heat flux data, which 
will be used in comparison later. From Fig. 3, it is clearly 
shown that the rewetting temperature profiles are transient in 
nature, even in the coordinate frame moving with the rewetting 
front. 
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The successful prediction of the rewetting curve by the ap- 
proximate analytical solution, Eq. (52), for the grooved (or 
smooth) plate with uniform heating is illustrated in Fig. 4. It 
is shown that the approximate closed-form solution is in rea- 
sonably good agreement with the numerical solution. This may 
appear to be somewhat of a surpise in view of the seemingly 
inconsistent approximation of treating P (or t/ r ) as constant. 
As mentioned above, such an approximation was necessarily 
made in the mathematical manipulation to achieve an ap- 
proximate closed-form solution. It is equivalent to neglecting 
higher order terms attributed to the transient components of 
P. Figure 4 shows that such an approximation does yield good 
results as is expected because all rewetting velocities tend to 
level off quickly. In fact, beyond t- 5.208 (or /= 1.4 s), the 
difference between the closed-form and the numerical solutions 
is also indiscernible. 

Figure 4 shows that the predicted rewetting velocity decreases 
with time. This is physically explainable because the plate tem- 
perature and, therefore, the thermal capacity to be removed 
ahead of the wet front increase with time due to heating beneath 
the plate. Consequently, the rewetting velocity slows down as 
it proceeds forward. In the case of quenching a hot plate 
without subjecting to a heating condition, as studied by pre- 
vious investigators, the rewetting velocity is constant. Thus, 
the present study reveals a fundamental difference between the 
quenching of a hot plate with heating and that without heating. 
In the former, the wetting speed is time dependent, while in 
the latter, it is simply time invariant. 

In a recent study, Grimley et al. (1988) conducted experi- 
ments to investigate the enhancement of convective boiling 
heat transfer by grooves on a plate heated from beneath. Un- 
fortunately no rewetting velocity data were reported that could 
otherwise be useful to check the validity of the solutions dis- 
cussed above. However, they did report interesting data on 
the maximum heat flux (the critical heat flux, CHF) that the 
heater could supply to the plate without causing dryout of the 
flowing liquid film. Since the dryout and the subsequent rewet 
of a heated surface are an integrated problem in heat pipe 
applications, it is desirable to be able to explain or predict the 
heat flux condition that leads to dryout of the plate. To achieve 
this objective, a simple method is presented next and com- 
parisons will be made with the reported data. 


Prediction of the Maximum Dryout Heat Flux 

Under consideration is a smooth or a grooved plate initially 
covered by a thin flowing liquid film at a temperature T s . The 
plate is subjected to a uniform heating. It is of interest to 
predict the maximum heat flux that triggers the dryout of the 
film. 

A simple method based on the above rewetting concept is 
now extended to provide a means of estimating this maximum 
dryout heat flux. For a flowing liquid Him over a plate heated 
by a heat flux that exceeds the maximum heat removal ca- 
pability by convection and boiling, the liquid will cease to 
advance and begin to recede. Thus dryout wilt occur. On this 
physical premise, the maximum rate of the heat removal is 
given by 

= (W) 


Following the work of Howard et al. (1975), the average con- 
vective boiling heat transfer coefficient of the liquid Him is 
estimated from 


-r 1 f 7 ”* Qt j T 

k To.,-T. ) r . T-T f dT 


(65) 


in which T, is the plate temperature at the onset of boiling, 
and T 0t , is the smooth plate Leidenfrost temperature of the 
rewetting front and is approximated by the plate temperature 
at the CHF location. Q b is the boiling curve of the liquid film 


Tabid 1 Comparison between predicted and measured (Grimley et 
al.,1988, Figs. 8 and 9) maximum dryout beat flux 
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over the plate. As an approximation, a form Qb = <*\(T- T,) n , 
which fits the boiling curve, can be used. In the event that the 
boiling curve of the flowing film is unavailable, the pool boiling 
curve could be used as the first-order approximation (Howard 
et al., 1975). 

In the case of the grooved plate, the same analogy developed 
above is adopted here. The maximum dryout heat flux is es- 
timated from 

q m3X , g = N.h.(T 0 , g -T s ) ( 66 ) 

where T 0 , g is the Leidenfrost temperature of the grooved plate, 
also approximated by its CHF temperature. Based on the ex- 
perimental data of Grimley et al. (1988) for a fluorocarbon 
(FC-72) liquid film falling over heated smooth and grooved 
plates, T 0t g was found to be slightly lower than T 0 ,,. Thus, if 
T 0 g of the grooved plate is unavailable due to the lack of 
boiling curve for the grooved plate, one may attempt to ap- 
proximate Tq, g from the smooth surface data To, Then the 
predicted maximum heat flux may be slightly overestimated, 
namely, 

q m >«. g <W-h-(T Q , s -T s ) 

As an application to show the feasibility of the above method , 
the experimental conditions and geometries of Grimley et al. 
(1988) are used. In their experiments, the liquid film completely 
covered the smooth and grooved plates, 6 « 0.5 mm, f 2 = 0.2 
mm, and *3 = 0.4 mm, such that N= 2.25 (from Eq. (58)). 
(7j — t s ) and (T 0 ~ T s ) are taken from their boiling curves. The 
average boiling heat transfer coefficient is calculated from their 
smooth surface boiling curve and Eq. (65). The two correlation 
constants a x and n are determined by arbitrarily collocating 
the boiling curve of the smooth surface plate at two locations, 
Tj and T 0 , where the boiling heat fluxes are designated by q, 
and <7 C hf» respectively. Table 1 shows the comparison between 
the predicted maximum dryout heat flux, q mzx , using Eqs (64) 
and (65), and the reported dryout data <7 C hf for both types of 
plates. The agreement is satisfactory particularly in view of 
the simplicity of the method proposed for the grooved plate. 
The same agreement is shown in Table 2 when the mean inlet 
velocity of the falling film is increased from 0.5 m/s to 1.0 
m/s. 

Conclusions 

The rewetting process of a smooth surface plate subject to 
a uniform heating has been investigated. A proper governing 
transient heat conduction with a convective boiling condition 
has been presented and solved to yield an approximate closed- 
form solution for the plate temperature profiles in the wet and 
dry regions of the plate. From the temperature profiles, an 
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approximate closed-form solution for the rewetting velocity 
over the heated plate has been obtained. Numerical solutions 
have also been presented to check the validity of the closed- 
form solution. The closed-form rewetting velocity was found 
to be in good agreement with that of the numerical solution. 
It is shown that in a limiting condition the present rewetting 
velocity solution reduces correctly to the existing solution for 
the rewetting of a hot, isothermal plate without heating. How- 
ever, contrary to the case without heating, the rewetting process 
on the plate with uniform heating is found to be transient (time 
variant) even on the coordinate frame moving with the re- 
wetting front. The rewetting velocity is found to be much faster 
initially and then levels off later. 

A method to address the rewetting process of the grooved 
plate based on the smooth plate rewetting model has been 
developed. It is shown that, by properly defined scalings, the 
solution for the smooth plate can be made to be applicable 
for the grooved plate. 

Finally, the dryout of a liquid film over a heat plate has 
been investigated. A simple method has been proposed to pre- 
dict the dryout critical heat flux of the smooth and grooved 
plates. The results of the prediction were compared and found 
to be in reasonable agreement with the existing experimental 
data. 
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APPENDIX A 
Coordinate Transformation 

The heat conduction equation within a smooth plate as shown 
in Fig. 1 is 


pCf 


(*t\ _r(£L. 

\dt) x , \dx ' 2 


£T 

a/ 


(Al) 


For a coordinate system x-o-y moving with the rewetting front, 


x =x + 


( U r (t)dt 

T(x\ t) = T(x{x\ t ), t) 


(A 2) 
(A3) 


Then, noting that 


be written as 

[©.-“©J-(S*S) “ 


Eq. (Al) can be written as 
pC, 


For a thin plate subjected to heating from below, it further 
reduces to 


= J^dri +u d_T 

l btl pC p bx 2 ' bx pC p 5, 


h a 

(T- T s ) + 


pC p S\ 


or in dimensionless form 

be b 2 e 


be 


T = iri +p ir-v 9+A 

07 or) or) 


(A 6) 


where the dimensionless variables are parameters are given by 
Eq. (2). 

For the purposes of comparison and numerical computation, 
the above is written in untransformed coordinates (x\ /) as 


where 


36 d 2 0 


7j' =x’/S] 


(A7) 

(A8) 
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